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Abstract 

In the hadronic decay mode of a pair of W bosons, e +e > W + W > ^1(12^3^45 

QCD interference effects can mix up the two colour singlets qiq2 and q3q4, i.e. produce 
hadrons that cannot be uniquely assigned to either of W + and W — . We show that 
interference is negligible for energetic perturbative gluon emission, and develop models to 
help us to estimate the non-perturbative effects. The total contribution to the systematic 
error on the W mass reconstruction may be as large as 40 MeV. 
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This letter is concerned with QCD interference effects that may occur when two unsta- 
ble particles decay and hadronize close to each other in space and time. W pair production 
is here of particular interest because of its practical importance and relative simplicity, 
but many other processes are also affected. (For further details and references see Ref. 

iO 

QCD interference effects between the W + and W — decays undermine the traditional 

meaning of a W mass in the process e -|-e + W > qiq2<l3Q4- Specifically, it 

is not even in principle possible to subdivide the hadronic final state into two groups of 
particles, one of which is produced by the qiq2 system of the W + decay and the other 
by the q3q4 system of the W — decay: some particles originate from the joint action of 
the two systems. Since a determination of the W mass is one of the main objectives of 
LEPlarge the ambiguities can be. A statistical error of 55the precision of the theoretical 
predictions should ideally match or exceed this experimental accuracy. 

A complete description of QCD interference effects is not possible since non- 
perturbative QCD is not well understood. The concept of colour reconnec- 
tion/rearrangement is therefore useful to quantify effects (at least in a first approxima- 
tion). In a reconnection two original colour singlets (such as qiq2 and q3q4) are transmuted 
into two new ones (such as qiq4 and q3q2)- Subsequently each singlet system is assumed 
to hadronize independently according to the standard algorithms, which have been so suc- 
cessful in describing e.g. Z decays. Depending on whether a reconnection has occurred 
or not, the hadronic final state is then going to be somewhat different. 

The colour reconnection effects were first studied by Gustafson, Pettersson and Zerwas 
0, but their results were mainly qualitative and were not targeted on what might actually 
be expected at LEPexample of the so-called instantaneous reconnection scenario, where 
the alternative colour singlets are immediately formed and allowed to radiate perturbative 
gluons. 

In order to understand which QCD interference effects can occur in hadronic W +W — 
decays, it is useful to examine the space-time picture of the process. Consider a typical 
cm. energy of 170 GeV, a W mass m-w = 80 GeV and width Fw = 2.08 GeV. The 
averaged (over the W mass distribution) proper lifetime for a W is (r) ~ (2/3)/i/Fw ~ 
0.06This gives a mean separation of the two decay vertices of 0.04A gluon with an energy 
uj ^ Fw therefore has a wavelength much smaller than the separation between the W + 
and W — decay vertices, and is emitted almost incoherently either by the qiq2 system 
or by the q3q4 one [Q. Only fairly soft gluons, u < Fw, feel the joint action of all four 
quark colour charges. On the other hand, the typical distance scale of hadronization is 
about Ithan the decay vertex separation. Therefore the hadronization phase may contain 
significant interference effects. 

In the following, we will first discuss perturbative effects and subsequently non- 
perturbative ones. (For a discussion of a possible interplay between the two stages see 
Ref. i.) 

Until today, perturbative QCD has mainly been applied to systems of primary partons 
produced almost simultaneously. The radiation accompanying such a system can be 
represented as a superposition of gauge-invariant terms, in which each external quark line 
is uniquely connected to an external antiquark line of the same colour. The system is thus 
decomposed into a set of colourless qq antennae/dipoles |^. One of the simplest examples 
is the celebrated qqg system, which (to leading order in 1/Nc 2, where Nq = 3 is the 
number of colours) is well approximated by the incoherent sum of two separate antennae, 
qg and gq. These dipoles radiate gluons, which within the perturbative scenario are the 
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principal sources of multiple hadroproduction. 

Neglecting interferences, the e +e > W + W > qiq2Q3Q^4 ^^^1 state can be 

subdivided into two separate dipoles, qiq2 and q3q4. Each dipole may radiate gluons 
from a maximum scale downwards. Within the perturbative approach, colour trans- 
mutations can result only from the interferences between gluons (virtual as well as real) 
radiated in the W + and W — decays. A colour reconnection then corresponds to radia- 
tion, e.g. from the dipoles qiq4 and q3q2. The emission of a single primary gluon cannot 
give interference effects, by colour conservation, so interference terms only enter in second 
order in as- 

The general structure of the results is well illustrated by the interference between the 
graph where a gluon with momentum ki (/C2) is emitted off the qiq2 (q3q4) dipole and the 
same graph with ki and k2 interchanged: 

(To uJi co'2 V 4vr 2 / Nc 2 — 1 

where Cp = {Nc 2 — l)/{2Nc) = 4/3. We proceed to comment on the non-trivial factors 
in this expression. 

The interference is suppressed by l/{Nc 2 — 1) = 1/8 as compared to the total rate of 
double primary gluon emissions. This is a result of the ratio of the corresponding colour 
traces. 

The so-called profile function xu @) controls decay-decay interferences. It quantifies 
the overlap of the W propagators in the interfering Feynman diagrams. Near the W +W — 
pair threshold, X12 simplifies to 

rw 2 

^'''^rw2 + (cui-cu2)2 ^ ' 

Other interferences (real or virtual) are described by somewhat different expressions, e.g. 
with ui — UJ2 —>■ Ui + UJ2, but have the same general properties. The profile functions 
cut down the phase space available for gluon emissions with u > Fw by the alterna- 
tive quark pairs. (We can neglect the contribution from kinematical configurations with 
uji,uj2 ^ Fw, — 1^21 < Fw since the corresponding phase-space volume is small.) The 
possibility for the reconnected systems to develop QCD cascades is thus reduced, i.e. the 
dipoles are almost sterile. 

The radiation pattern H{k) is given by 

H{k) = qiq4 + qaqs - qlqs - ^^4 (3) 

where the radiation antennae are [^] 

fj = ■ (4) 

{pi ■ k){pj ■ k) 

In addition to the two dipoles qiq4 and q3q2, which may be interpreted in terms of 
reconnected colour singlets, one finds two other terms, q^qs and q2q4, which come in with 
a negative sign. The signs represent the attractive and repulsive forces between quarks 
and antiquarks [|], 0]. 

It should be emphasized that, analogously to other colour-suppressed interference 
phenomena, rearrangement can be viewed only on a completely inclusive basis, when all 
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the antennae are simultaneously active in the particle production. The very fact that 
the reconnection pieces are not positive-definite refiects their wave interference nature. 
Therefore the effects of reconnected almost sterile cascades should appear on top of a 
dominant background generated by the ordinary-looking no-reconnection dipoles qiq2 
and q3q4. 

Summing up the above discussion, it can be concluded that perturbative colour recon- 
nection phenomena are suppressed, firstly because of the overall factor ^/{Nc 2—1), and 
secondly because the rearranged dipoles can only radiate gluons with energies u < Fw 
Only a few low-energy particles should therefore be affected, AA^ recon/A^ no — recon < 
C(10 -2). 

We now turn to the possibility of reconnection occurring as a part of the non- 
perturbative hadronization phase. Since hadronization is not understood from first princi- 
ples, this requires model building rather than exact calculations. We will use the standard 
Lund string fragmentation model |^ as a starting point, but have to extend it consid- 
erably. The string is here to be viewed as a Lorentz covariant representation of a linear 
confinement field. 

The string description is entirely probabilistic, i.e. any negative-sign interference ef- 
fects are absent. This means that the original colour singlets qiq2 and q3q4 may transmute 
to new singlets qiq4 and q3q2, but that any effects e.g. of qiqa or q2q4 dipoles are absent. 
In this respect, the non-perturbative discussion is more limited in outlook than the per- 
turbative one above. However, note that dipoles such as qiqa do not correspond to colour 
singlets, and can therefore not survive in the long-distance limit of the theory, i.e. they 
have to disappear in the hadronization phase. 

The imagined time sequence is the following. The W + and W — fly apart from their 
common production vertex and decay at some distance. Around each of these decay ver- 
tices, a perturbative parton shower evolves from an original qq pair. The typical distance 
that a virtual parton (of mass m ~ lOhadronic flnal state) travels before branching is 
comparable with the average W +W — separation, but shorter than the fragmentation 
time. Each W can therefore effectively be viewed as instantaneously decaying into a string 
spanned between the partons, from a quark end via a number of intermediate gluons to 
the antiquark end. The strings expand, both transversely and longitudinally, at a speed 
limited by that of light. They eventually fragment into hadrons and disappear. Before 
that time, however, the string from the W + and the one from the W — may overlap. If 
so, there is some probability for a colour reconnection to occur in the overlap region. The 
fragmentation process is then modifled. 

The Lund string model does not constrain the nature of the string fully. At one 
extreme, the string may be viewed as an elongated bag, i.e. as a flux tube without any 
pronounced internal structure. At the other extreme, the string contains a very thin core, 
a vortex line, which carries all the topological information, while the energy is distributed 
over a larger surrounding region. The latter alternative is the chromoelectric analogue 
to the magnetic flux lines in a type II superconductor, whereas the former one is more 
akin to the structure of a type I superconductor. We use them as starting points for two 
contrasting approaches, with nomenclature inspired by the superconductor analogy. 

In scenariospace-time volume over which the W + and W — strings overlap, with 
saturation at unit probability. This probability is calculated as follows. In the rest frame 
of a string piece expanding along the ±z direction, the colour fleld strength is assumed 
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to be given by 

Q{^, t) = exp {-{x 2 + y 2)/2rhad 2} 9{t - |x|) exp {-{t 2 - z 2)/rfr,g 2} (5) 

The first factor gives a Gaussian fall-off in the transverse directions, with a string width 
^had ~ 0.5of typical hadronic dimensions. The time retardation factor 6(t — |x|) ensures 
that information on the decay of the W spreads outwards with the speed of light. The last 
factor gives the probability that the string has not yet fragmented at a given proper time 
along the string axis, with r^ag ~ l.Sfrom the W + decay, this field strength has to be 
appropriately rotated, boosted and displaced to the W + decay vertex. In addition, since 
the W + string can be made up of many pieces, the string field strength fimax + (x, t) is 
defined as the maximum of all the contributing +'s in the given point. The probability 
for a reconnection to occur is now given by 

Brecon = 1 - exp ^-A;i j d 3x dt r^max + (x, t) fi^ax - (x, t)^ (6) 

where ki is a free parameter. If a reconnection occurs, the space-time point for this 
reconnection is selected according to the differential probability Qmax +(x, t) fimax — (x, t). 
This defines the string pieces involved and the new colour singlets. 

In scenario II it is assumed that reconnections can only take place when the core 
regions of two string pieces cross each other. This means that the transverse extent of 
strings can be neglected, which leads to considerable simplifications compared with the 
previous scenario. The position of a string piece at time t is described by a one-parameter 
set x(t, a), where < a < 1 is used to denote the position along the string. To find 
whether two string pieces i and j from the W + and W — decays cross, it is sufficient to 
solve the equation system Xj + (t, a +) = x^ — (t, a — ) and to check that this (unique) 
solution is in the physically allowed domain. Further, it is required that neither string 
piece has had time to fragment, which gives two extra suppression factors of the form 
exp{— r 2/rfi.ag 2}, with r the proper lifetime of each string piece at the point of crossing, 
i.e. as in scenariostring crossings, only the one that occurs first is retained. 

Both scenarios are implemented in a detailed simulation of the full process of W ± 
production and decay, parton shower evolution and hadronization 0. It is therefore 
possible to assess any experimental consequences for an ideal detector. 

The reconnection probability is predicted in scenariothe possibility to vary the baseline 
model in a few respects. Scenariofci. We have chosen ki to give an average Prccon ~ 0.35 
at 170 GeV, as is predicted in scenario 

The resulting cm. energy dependence of Precon is very slow: between 150 and 200a 
factor of 2. Here it is useful to remember that the W ± are never produced at rest with 
respect to each other: the naive Breit-Wigner mass distributions are distorted by phase- 
space effects, which favour lower W masses. For 150-200in the range 22-60 GeV, rather 
than in the range 0-60It is largely this momentum that indicates how fast the two W 
systems are flying apart, and therefore how much they overlap in the middle of the event. 
Also the energy variation in the perturbative description is very small. If we want to call 
colour reconnection a threshold effect, we have to acknowledge that the threshold region 
is very extended. 

Comparing the scenarios I and II above with the no-reconnection scenario, it turns out 
that reconnection effects are very small. The change in the average charged multiplicity 
is at the level of a per cent or less, and similar statements hold for rapidity distributions. 
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thrust distributions, and so on. This is below the experimental precision we may expect, 
and so may well go unobserved. One would like to introduce more clever measures, which 
are especially sensitive to the interesting features, but so far we have had little success. 

Ultimately, the hope would be to distinguish scenariosinto the nature of the confine- 
ment mechanism. In principle, there are such differences. For instance, the reconnection 
probability is much more sensitive to the event topology in scenarioof having two string 
cores cross is more selective than that of having two broad flux tubes overlap. 

We now come to the single most critical observable for LEPphysics, namely the W 
mass. Experimentally, depends in a non-trivial fashion on all particle momenta of an 
event. Errors in the W mass determination come from a number of sources 0, which we 
do not intend to address here. Therefore we only study the extent to which the average 
reconstructed W mass is shifted when reconnection effects are added, but everything else 
is kept the same. Even so, results do depend on the reconstruction algorithm used. We 
have tried a few different ones, which however all are based on the same philosophy: a jet 
flnder is used to deflne at least four jets, events with two very nearby jets or with more 
than four jets are rejected, the remaining jets are paired to deflne the two W's, and the 
average W mass of the event is calculated. Events where this number agrees to better 
than 10 GeV with the input average mass are used to calculate the systematic mass shift. 

In scenario I this shift is consistent with being zero, within the 10(160,000 events 
per scenario). Scenarioshift, of about — 30of the basic scheme. A simpler model, where 
reconnections are always assumed to occur at the centre of the event, instead gives a 
positive mass shift: about +30Precon ~ 0.35. We are therefore forced to conclude that 
not even the sign of the effect can be taken for granted, but that a real uncertainty of 
±30from our ignorance of non-perturbative reconnection effects. 

To examine the perturbative rearrangement effect, we have used a scenario where the 
original qiqg and q3q4 dipoles are instantaneously reconnected to qiq4 and q3q2 ones, and 
these are allowed to radiate gluons with an upper cut-off given by the respective dipole 
invariant mass . This gives a mass shift by about -|-500argued that real effects would be 
suppressed by at least a factor of 10 —2 compared to this, and thus assign a 5this source. 
Finally, the possibility of an interplay between the perturbative and non-perturbative 
phases must be kept in mind. We have no way of modelling it, but believe it will not be 
much larger than the perturbative contribution, and thus assign a further Sindependent, 
the numbers are added linearly to get an estimated total uncertainty of 40 

In view of the aimed-for precision, 40non- negligible. However, remember that as a 
fraction of the W mass itself it is a half a per mille error. Reconnection effects are 
therefore smaller in the W mass than in many other observables, such as the charged 
multiplicity. 

Clearly, it is important to study how sensitive experimental mass reconstruction algo- 
rithms are, and not just rely on the numbers of this paper. We believe that the uncertainty 
can be reduced by a suitable tuning of the algorithms, e.g. with respect to the impor- 
tance given to low-momentum particles, and with respect to the statistical treatment of 
the wings of the W mass distribution. The 40to be assigned to any algorithm that has 
not been properly evaluated. 

In summary, we have developed the flrst detailed model of QCD rearrangement effects 
in the decay of two heavy colourless objects into quarks and gluons. Beyond the immediate 



use for LEPother processes, e.g. top quark production and decay H, 0, 10 
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